A non-standard symmetry-preserving method to compute bounded solutions of a generalized Newell-Whitehead-Segel equation

In this work, we propose a finite-difference scheme to approximate the solutions of a generalization of the classical, one-dimensional, Newell-Whitehead-Segel equation from fluid mechanics, which is an equation for which the existence of bounded solutions is a well-known fact. The numerical method preserves the skew-symmetry of the problem of interest, and it is a non-standard technique which consistently approximates the solutions of the equation under investigation, with a consistency of the first order in time and of the second order in space. We prove that, under relatively flexible conditions on the computational parameters of the method, our technique yields bounded numerical approximations for every set of bounded initial estimates. Some simulations are provided in order to verify the validity of our analytical results. In turn, the validity of the computational constraints under which the method guarantees the preservation of the boundedness of the approximations, is successfully tested by means of computational experiments in some particular instances.     

On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows

We present a new unconditionally positivity‐preserving (PP) implicit time integration method for the DG scheme applied to shallow water flows. This novel time discretization enhances the currently used PP DG schemes, because in the majority of previous work, explicit time stepping is implemented to deal with wetting and drying. However, for explicit time integration, linear stability requires very small time steps. Especially for locally refined grids, the stiff system resulting from space discretization makes implicit or partially implicit time stepping absolutely necessary. As implicit schemes require a lot of computational time solving large systems of nonlinear equations, a much larger time step is necessary to beat explicit time stepping in terms of CPU time. Unfortunately, the current PP implicit schemes are subject to time step restrictions due to a so‐called strong stability preserving constraint. In this work, we hence give a novel approach to positivity preservation including its theoretical background. The new technique is based on the so‐called Patankar trick and guarantees non‐negativity of the water height for any time step size while still preserving conservativity. In the DG context, we prove consistency of the discretization as well as a truncation error of the third order away from the wet–dry transition. Because of the proposed modification, the implicit scheme can take full advantage of larger time steps and is able to beat explicit time stepping in terms of CPU time. The performance and accuracy of this new method are studied for several classical test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

On the fluidization of Petri nets and marking homothecy

The analysis of Discrete Event Dynamic Systems suffers from the well known state explosion problem. A classical technique to overcome it is to relax the behavior by partially removing the integrality constraints and thus to deal with hybrid or continuous systems. In the Petri nets framework, continuous net systems (technically hybrid systems) are the result of removing the integrality constraint in the firing of transitions. This relaxation may highly reduce the complexity of analysis techniques but may not preserve important properties of the original system. This paper deals with the basic operation of fluidization. More precisely, it aims at establishing conditions that a discrete system must satisfy so that a given property is preserved by the continuous relaxation. These conditions will be mainly based on the marking homothetic behavior of the system. The focus will be on logical properties as boundedness, B-fairness, deadlock-freeness, liveness and reversibility. Furthermore, testing homothetic monotonicity of some properties in the discrete systems is also studied, as well as techniques to improve the quality of the fluid relaxation by removing spurious solutions.     

模糊互补判断矩阵的一个通用排序公式及其保序性研究

基于模糊互补判断矩阵的一致性转化,利用行和归一化方法,得到了模糊互补判断矩阵的一个排序公式,指出目前基于模糊加性一致的排序公式大多都是所提方法的特例,并给出了在增加一个或一组新元素时的保序性条件.     

Development of level set method with good area preservation to predict interface in two‐phase flows

A two‐step conservative level set method is proposed in this study to simulate the gas/water two‐phase flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the coupled compact scheme. For accurately predicting the modified level set function, the dispersion‐relation‐preserving advection scheme is developed to preserve the theoretical dispersion relation for the first‐order derivative terms shown in the pure advection equation cast in conservative form. For the purpose of retaining its long‐time accurate Casimir functionals and Hamiltonian in the transport equation for the level set function, the time derivative term is discretized by the sixth‐order accurate symplectic Runge–Kutta scheme. To resolve contact discontinuity oscillations near interface, nonlinear compression flux term and artificial damping term are properly added to the second‐step equation of the modified level set method. For the verification of the proposed dispersion‐relation‐preserving scheme applied in non‐staggered grids for solving the incompressible flow equations, three benchmark problems have been chosen in this study. The conservative level set method with area‐preserving property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam‐break, Rayleigh–Taylor instability, bubble rising in water, and droplet falling in water problems. Good agreements with the referenced solutions are demonstrated in all the investigated problems. Copyright © 2010 John Wiley & Sons, Ltd.     

A multidimensional HLL‐Riemann solver for non‐linear hyperbolic systems

This article presents a numerical model that enables to solve on unstructured triangular meshes and with a high order of accuracy, Riemann problems that appear when solving hyperbolic systems. For this purpose, we use a MUSCL‐like procedure in a ‘cell‐vertex’ finite‐volume framework. In the first part of this procedure, we devise a four‐state bi‐dimensional HLL solver (HLL‐2D). This solver is based upon the Riemann problem generated at the barycenter of a triangular cell, from the surrounding cell‐averages. A new three‐wave model makes it possible to solve this problem, approximately. A first‐order version of the bi‐dimensional Riemann solver is then generated for discretizing the full compressible Euler equations. In the second part of the MUSCL procedure, we develop a polynomial reconstruction that uses all the surrounding numerical data of a given point, to give at best third‐order accuracy. The resulting over determined system is solved by using a least‐square methodology. To enforce monotonicity conditions into the polynomial interpolation, we use and adapt the monotonicity‐preserving limiter, initially devised by Barth (AIAA Paper 90‐0013, 1990). Numerical tests and comparisons with competing numerical methods enable to identify the salient features of the whole model. Copyright © 2010 John Wiley & Sons, Ltd.     

变截面主流加速对超音速气膜冷却的影响

本文对扩张通道和直通道中的超音速气膜冷却进行了数值模拟,分析了主流加速和不加速两种情况下超音速气膜冷却的差异.计算结果表明,在主流加速的情况下,超音速气膜冷却的效果要好于相同进口条件下主流不加速的情况,并且这种差异越到下游越明显.同时,在本文计算的工况范围中,无论主流加速还是不加速,增大主流进口马赫数或者减少冷却流进口马赫数都使气膜冷却效率下降,而壁面静压则由主流气体决定.     

General Theory of Global Smoothness Preservation by Singular Integrals, Univariate Case

In this article it is established that the well-known singular integrals of Picard, Poisson–Cauchy, Gauss–Weierstrass, and their Jackson-type generalizations fulfill the global smoothness preservation property, i.e., they ripple less than the function they are applied on, that is, producing a nice and fit approximation to the unit. The related results are established over various spaces of functions and the associated inequalities involve different types of corresponding moduli of smoothness. Several times these inequalities are proved to be sharp, namely, they are attained.